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The relation between the baryon number in QCD at nonzero chemical potential and the 
spectral density of the baryon number Dirac operator, jo(D + m), is examined. We show that 
extreme oscillations of the spectral density, caused by the QCD sign problem, are essential 
for the formation of the average baryon number when [i > m w /2. We compute the oscillating 
region of the spectral density using chiral perturbation theory. The extreme oscillations have 
a microscopic period and are resolved using random matrix theory. 
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I. INTRODUCTION 

Lattice QCD [1] is our best nonperturbative tool to study the phase structure of strongly interacting 
matter as a function of temperature. At low temperature the formation of the chiral condensate, 
which signals the spontaneous breaking of chiral symmetry, can be read of from the density of 
eigenvalues of the lattice Dirac operator at the origin. This useful link between the chiral condensate 
and the eigenvalues of the Dirac operator is known as the Banks-Casher relation ^]. 

At nonzero chemical potential, /z, the Monte Carlo method, which forms the basis of Lattice 
QCD simulations, is hampered by the sign problem: The fermion determinant of the Dirac operator, 
and consequently the measure on which we would like to perform Monte Carlo sampling, takes 
complex values, see for example the reviews [Si. Also the Banks-Casher relation 2fl breaks down 

nn 

at nonzero chemical potential. The new relation |5(, which replaces the Banks-Casher relation, 
explicitly shows the intimate connection between the sign problem and the spontaneous breaking 
of chiral symmetry: it is a region with extreme oscillations of the eigenvalue density of the Dirac 
operator [g,Q], caused by the sign problem, which leads to the discontinuity of the chiral condensate 
in the chiral limit at nonzero chemical potential. The oscillations have a period of order the inverse 
volume and an amplitude which is exponentially large in the volume. The region within which the 
oscillations take part can be computed using the mean field approximation in chiral perturbation 



theory 



but as the period of the oscillations are on the microscopic scale, set by the inverse 



volume, one needs the exact microscopic eigenvalue density jf], Q] in order to resolve them. For a 
review see e.g. {{J. 

The new relation 0,2], which replaces the Banks-Casher relation, not only give a direct insight 
in the way the chiral condensate is formed in unquenched QCD at nonzero chemical potential. It 
also solves the long standing problem of how the eigenvalues of the Dirac operator, D, can have 
a strong dependence on the chemical potential while at low temperature the partition function is 
independent of the chemical potential when [i is less than a third of the nuclon mass (in [10|] this 
was coined the Silver Blaze problem) . The relation between a region of the eigenvalue density with 
extreme oscillations and the spontaneous breaking of chiral symmetry has also been established in 
ldQCD 



111 ] as well as in two color QCD with non degenerate quark masses [121 ] . 
Here we show that the mechanism behind the new relation [J, |5] is not restricted to the chiral 
condensate and the spectral density of D: It is the exact same mechanism that links the baryon 
number density to the spectral density of baryon number Dirac operator, jo(D + m). 
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The average baryon number is given by 

Vn B (jj) = ^-logZ( M ) = (Tr 1 ), (1) 
a/i \ L>o(m) + fi/ 

where the baryon number Dirac operator is, Do(m) = 7o(-D + m). The average trace in ([T]) can be 
expressed as the integral 

Vn B (fx)= [ -^Lp(A,A*), (2) 

where A are the eigenvalues of the baryon number Dirac operator and p is the spectral density of 
the baryon number Dirac operator. 

We will compute the eigenvalue density of Do(m) and demonstrate that it is again a strongly 
oscillating region of the unquenched eigenvalue density which ensures the correct physical behavior 
of the average baryon density. In order to establish this result we will investigate the structure of 
the spectral density of the baryon number Dirac spectrum using chiral perturbation theory (%PT) 
and random matrix theory (RMT). 

In quenched simulations a strongly oscillating region of the eigenvalue density is not possible, 
since the quenched theory is free from the sign problem. In such simulations the baryon number 
density therefore becomes nonzero when the chemical potential reaches half the pion mass rather 
than at a third of the baryon mass as expected for unquenched QCD. This early onset of the baryon 
density in the quenched theory has been observed on the lattice Q and in a U(l) model 
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It was also understood analytically 16] within a random matrix model (RMM). 

Before we turn to the actual computation of the eigenvalue density of the baryon number Dirac 
operator we give, in section HH a simple example which illustrates the mechanism involving a 
strongly oscillating density. Then in section IIIII we will use a mean field approach to describe 
the phase diagram of the spectrum. Finally, in section IIV| we will use a RMM to show that the 
extreme oscillations of the microscopic eigenvalue density lead to the expected behavior of the 
average baryon number density. 

II. AN EXAMPLE WITH OSCILLATING REGIONS 

In this section we will give a simple example to illustrate the point that oscillations with a 
microscopic wavelength can remove the unphysical early onset of the average baryon number at 
fi = m n /2 observed in the quenched theory. 

First, let us consider the nonoscillating eigenvalue density corresponding to the quenched case 

m*) = \f%v(i + - m n /2), (3) 
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where x and y are the real and imaginary part of the eigenvalue A, see figure [TJ This spectral 
density results, through |2|, in a baryon density, which is zero for /i < m^/2 and nonzero for 

H > m vr /2, 

= F$V J dxsign(x + /i)(l + ^)#(|x|-mj2) (4) 

= 2^(/x-^)^(H-m 7r /2). 

This is like the observed quenched baryon density, see figure [2j Let us now illustrate the main 
point: The unquenched average baryon number is zero also when /i is in the range between half 
the pion mass and a third of the nucleon mass because of two strongly oscillating regions. In order 
to see how this works we extend the example as 

pVf(x, y, M , = p% x (x, y, m^)(l - e F ' v ^>^)), (5) 

with 

U(x,y,[i,m) = -{y + isign(ar)(|ar| + |^|)) 2 - ml - (4x 2 - ml)9{\x\ - m w /2), (6) 

and where V is the four volume. Note that the amplitude of the oscillations grows exponentially 
with V, and that the period of the oscillations (determined by e~ 2iysignx ^ x ^ tJ, ^ F ^ v ) is of order 
1/V. 

The contribution to the average baryon density from the oscillating part of the eigenvalue density 
can be evaluated by a saddle point integration, in the large V limit. The contour along the real y 
axis is deformed to go through the saddle point y = — i sign(x)(|x| + |/x|). If the deformation passes 
the pole at y = i(x + fi), then we add a contribution for integration around the pole. We obtain 

roo F%VU(x,y,fi,m) 

/ dy — — = 2*[0{ji)0{ r x)0{x + ii)-0{-i i )0(x)0(-x-ii)]. (7) 

J-oo X ~T ly ~T [I 

It follows that the average baryon number is zero also for fi > m^/2 

Vn?>™>M = J _^-^(x iyiM) m.) = 0. (8) 

What we have learned from this example is that bounded oscillating regions can remove the un- 
physical early onset of the average baryon number at [i = m n /2. Below we show, using xPT and 
RMT, that this mechanism is realized in unquenched QCD at nonzero [i. 
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III. THE BOUNDARIES OF THE OSCILLATING REGIONS 

In this section we will investigate the phase diagram of the baryon number Dirac spectrum by 
means of xPT at mean field level. Our aim is to calculate the boundaries of the oscillating regions. 
To do this we employ the replica trick (see e.g. [13]), where n replica pairs consisting of quarks and 
conjugate quarks has been introduced in the partition function. The replica partition function is 
given by the ensemble average of Nf quarks and n replica pairs, 

Z N f' n (\,\*,fj,,m) = (det n (D (m) + X)det n (D (m) - X*)dei N f (D (m) + /A (9) 

The spectral density of the baryon number Dirac operator is obtained from the replica partition 
function by (see [la ] for a detailed explanation of how the density results from the derivatives) 

p N f(X,X*,H,m) = lim —d\d\* log Z N f' n (X, A*, u, m). (10) 

n-5-0 n 

Note that the chemical potential of the replica quarks, A, corresponds to the eigenvalues of -Do( m )- 
We will examine the partition function Z N t' n (X, A*, //, m) to leading order in chiral perturbation 
theory. The chiral Lagrangian to leading order is given by 3] 

C cS = -^-TrD u UD u U^ - |TrM(C/ + C/t) ; (n) 

where X is the magnitude of the chiral condensate and F n is the pion decay constant. The chemical 
potential enters the QCD partition function as en external vector field, B v = 5 U $B. To respect 
local invariance, B v can only appear in the covariant derivative, see e.g. 
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D U U = d u U + i[U,B u \. (12) 

In the case where we have Nf quarks and n replica pairs all with mass m and with chemical 
potential [i and A, respectively, then we have that 

U G SU (2n + Nf), M = ml 2n+Nf and B = diag(AI n , -A*I n , fd N/ ). (13) 

A. The Quenched Theory 

Here we look at the quenched theory, where Nf = 0, and the partition function ([9]) is given by 
the quenched ensemble average of n replica pairs, 

Z n (X, X*,m) = <det n (D (m) + A) det n (D (m) - A*)), (14) 



FIG. 1: Left panel: The phase diagram of the quenched baryon number Dirac spectrum in the complex 
eigenvalue A-plane. The diagram has a strip of zero spectral density and width m v centered around the 
imaginary A axis. The normal phase is labeled by N and n refers to pion condensed phase. The eigenvalue 
density is nonzero in the pion condensed phase. Right panel: Plot of the quenched spectral density of the 
baryon number Dirac operator as a function of Re A (the density is independent of Im A). If m — 0, the pion 
is massless and the spectral density is represented by the dotted line. 

where m is the mass of the replica quarks and A is the complex chemical potential. The spectral 
density of the baryon number Dirac operator is obtained as 



In the mean field limit the free energy has a trivial linear dependence on n. For this reason we can 
restrict ourselves to n = 1 in the following calculation. From ()14[) with n = 1 one sees that the 
chemical potential of replica quarks corresponds to isospin chemical potential for real A. At low 
temperature and ji > m n /2 there will be a condensate of quarks and conjugate quarks. One can 
interpret the quarks and conjugate quarks as "up" and a "down" quarks, and for this reason we 
will refer to the condensate of replica quarks as the pion condensate. 

The partition function Z n=1 (A, X*,m) depends on the chiral Lagrangian given in (llip and (I12D 
with U £ SU{2), M = diag(m, m) and B = diag(A, —A*). The mean field phase diagram, however, 
only depends on the static part of the Lagrangian 




n->0 n 



(15) 



C, 



•stat 




TtM(U + U ] ). 



(16) 



To find the mean field structure we use the ansatz 




(17) 
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FIG. 2: A plot of the quenched mean field baryon number density (solid line). If m = 0, then the pions are 
massless and baryon density is represented by the dotted line. 



which results in an effective Lagrangian 



(18) 



The minima of this Lagrangian are easily found to be at sina = and cos a = 2m£/i ?2 (A + A*) 2 . 
This yields the free energies 

■.N 



vc 



stat 



v ''-stat 



-2mT,V, 

F*V 2 2m 2 SV 2 

2 ( + ' F*V(\ + A*) 2 ' 



(19) 



for the normal and pion phase, respectively. The boundary between the two phases is given by 

2mE 



(Re A)" 



4F2 



2 



(20) 



where we used the Gell-Mann-Oakes-Renner relation for the Goldstone mass [23], to 2 = 2mS/i ?2 . 
Using (|15p we find that the spectral density is zero in the normal phase and given by 

12m 2 T, 2 V ~ / 3ml 



1 + 



in the pion phase. This structure is consistent with that obtained from a RMM in 



(21) 



24 



25(. Note 



that the phases of the replicated partition function are directly linked to the behavior of the 
eigenvalue density, hence we speak of the phase diagram of the baryon number Dirac spectrum. In 
terms of the microscopic parameters rh = m'EV and A = XF^^/V we have 



p*(x,rh) = F 2 V 1 + 



3rh 2 
4 ~¥ 



(22) 
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where x = Re A. A plot of the quenched phase diagram and the spectral density is given in figure[TJ 
We will now look at the quenched average baryon number, 

Vn^ Q {n,m) = [ d 2 A P ^ A ' A *' m) , (23) 

where A is the eigenvalues of Do(m), and the quenched mean field spectral density is given by 

p« F (A, A*,m) = F$v(l + \%)0{\x\ - mj2). (24) 

This is precisely the quenched density we used in the example and again integrating over the real 
and imaginary part of A leads to 




On figure [2] one sees that the quenched theory predicts a nonzero baryon density for \i > m^/2 as 
mentioned in the introduction. For m = the mass of the Goldstone bosons are zero. In this case 
the spectral density is a constant and the average baryon density is linear (see figure [2]). 



B. The Unquenched Theory 

We now consider the unquenched case. Just as in the queched theory the phases of the replicated 
partition function are directly linked to the behavior of the eigenvalue density of the baryon number 
Dirac opartor. In mean field the boundaries between the phases in the unquenched phase 
diagram of the baryon number Dirac spectrum are independent of the number of quark flavors and 
replicas. For this reason we will concentrate on the case n = 1 and Nf = 1. The replica partition 
function is given by 

Z N f =1 > n=1 {\, A*,//, m) = /det(A)(m) + A) det(£> (m) - A*) det(D Q (m) + n)\, (26) 

where m is the mass of both ordinary quarks and replicas, \x is the ordinary chemical potential 
and A is the replica chemical potential. Because of the similarity with the phase diagram of the 
(u,d,s) quark triplet j^fj, we will refer to the condensate of ordinary quarks with replicas as the 
kaon condensate. The masses of the Goldstone bosons (pion and kaon) are in this case equal 

ml = m\ = 2mJ:/Fl (27) 
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The static chiral Lagrangian to lowest order is the same as in equation (|16p but with M = 777JI3, 
B = diag(A, — A*, fj,) and of XJ € SU(3). In this case we make an ansatz for the Goldstone fields 
given by 

U = R 1 (a)R 3 ^), (28) 

where Rj(a) is a rotation by a about the j axis. This implies that the static Lagrangian becomes 
F 2 

r — —Ljl 

•'-stat — ^ 



2(A* + nf + (A + A*) (A* - A + 2/i) sin^ (3 ) sin 2 a + (A + A*) 2 sin 2 /3(1 + cos 2 a) 

-mS(cos/3 + cosa(l + COS/3)). (29) 



We will focus on the minima with /i > and x = Re A > 0, since the rest of the phase diagram is 
given by symmetry. For a = there are two minima: 

sin/3 = and cos (3 = 2mT./F 2 (\ + A*) 2 . (30) 

These minima give the free energy of the normal and pion condensed phases 

V£%* = -3mSV, 

stat 2 F%V(\ + A*) 2 v ~ ' 77 ' 

Note that the A and A* dependence at these minima are exactly as in the quenched case, cf. (|19p . 
When these minima are dominant the density is therefore identical to that found in the quenched 
case. 

The unquenched theory has a third minimum given by a = vr/2 and cos/3 = 2mS/F 2 (A* + /x) 2 , 
with implies the free energy 

VC stat - mW - F 2 V{fl + X*)2- ^ lb ) 

The phase boundaries occur where the real part of the free energy of two different phases are equal. 
The boundary between the normal and pion condensed phase is given by half the pion mass, just as 
in the quenched case. The formal expression of the normal-kaon and the pion-kaon boundaries are 
both given by a cubic equation in y 2 = Re A 2 with one real root and two complex conjugate roots. 
The structure of the kaon boundary can be seen on figure [3] and [H For /i < m n /2 the unquenched 
phase diagram is identical to the quenched case (see figured]). This is exactly what is expected, 
since the quenched theory trivially gives correct predictions for the average baryon density as long 
as \i < m n /2. For fi > the kaon phase appears in the unquenched spectral phase diagram, 
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Y Re A 

FIG. 3: Unquenched phase diagram in the Re A — \i plane. The labels N, it and K refer to the normal, pion 
and kaon condensed phases. The dashed line markes the boundary between the normal and kaon phase. 
To the right of the full line in the kaon phase the eigenvalue density is expected to be strongly oscillating. 
Note that the kaon condensed phase is introduced in the spectrum at /i = which is exactly where the 

unphysical behavior of the quenched theory occur. 

see Figure O We expect that the eigenvalue density is strongly oscillating within the kaon region 
and that these oscillations cure the unphysical behavior of the quenched baryon density, just as in 
the example in section |IlJ A diagram which shows the phase boundaries for fi > m n /2 of the kaon 
condensed phase are given in figure HI The details of the oscillations of the eigenvalue density are 
not resolved at the mean field level because the period of the oscillations is of order 1/V. Rather 
at mean field level one finds simply p = in the kaon phase. In section IIV Al we will compute the 
microscopic eigenvalue density using a RMM and show that the kaon condensed phase is indeed a 
strongly oscillating region, which exactly cancels the unphysical early onset of the average baryon 
density observed in the quenched case. 

C. The Massless Case 

In the limit of zero quark mass an additional saddelpoint is present. Here we analyze this case. 
In section IIV Al will compute the microscopic spectral density using a RMM of massless quarks, 
and the results of the present subsection will allow us to compare directly between %PT an d the 
RMM. 
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N 





FIG. 4: Left panel: Phase diagram for the spectrum of the baryon number Dirac operator with /i > 

The labels N, it and K refer to the normal, pion and kaon condensed phases. The density is expected to be 



strongly oscillating within the regions marked by K, which crosses the real axis at x = ±j(— M+V / i2 + 2m^) 
and x — ±/i. Right panel: Phase diagram for the spectrum of the baryon number Dirac operator with 
massless quarks (m = 0). The circular region contain strong oscillations (see section HV A[) and has radius 
\i. 

At m = the static part of the chiral Lagrangian (|16p is given by 
F 2 



r — —Ul 

*-\stat — 



2(A* + nf + (A + A*)(A* - A + 2/i) sm z (3 ) sin z a + (A + A*) 2 sin z (3(1 + cos z a) 



(32) 



which has four minima. The minima and corresponding free energies are 

Phase a (3 Free energy, VC (a, f3) 



N 
K 

k 



o 









2 







-2x 2 

2 " -\{x-iy + fj) 2 

f f i(_3x 2 + (y + iA) 2 ) 

The first three free energies can be found from the massive mean field results (|3ip by setting m = 0; 
the additional minimum of the massless theory exist since all minima trivially fulfill the constraints 
|cosa| < 1 and |cos/3| < \. 

The free energy of the normal phase is never a global minimum in the theory of massless 
quarks. The phase diagram of the quenched theory is therefore everywhere dominated by the 
condensate of the (massless) pions. Unquenching introduces condensation between replica and 
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physical flavors. The phase boundary between these phases and the pion condensed phase can 
be seen in the right hand panel of figure S] as the solid lines. The circular region is given by the 
boundary between pion and the K phase, and in section IIV Al we will show that this phase is 
dominated by strong oscillations, which exactly cancels the unphysical behavior of the quenched 
theory. In case of massless quarks the kaon condensed phase is a false minimum coming from the 
incomplete description of the fermionic replica trick. For a critical discussion of the replica method 



sec 



27 



281 ] . The problem becomes explicit within the supersymmetric technique (see 29|, |30] for 



an introduction) where one obtains the density from 



p N f (X, X* , fx,m) = lim d x d x *logZ N f +2 \ 2 (X,X*,X',\'*,fM,m), (33) 

A'— ¥\ 



where 



7 7V f +2|2 n x * w w* s / det(£> (m) + A)det(A)(m) -A*) N \ 

Zf (A, A ,A,A ^^) = ( v detpo(m) + AOdet(Z)o(m) _ A . ) det ' (D (m) + A }. (34) 

The presence of both fermionic and bosonic determinants leads to both ordinary and fermionic 
Goldstone modes in the effective xPT formulation. Integrating out the noncommuting variables 
gives a partially quenched partition function of the form 

Z N f +2 ^ 2 = J dfi F J dn B Pe- SF e~ SB , (35) 

where dfiF {dps) is the integration measure of the commuting variables related to the fermionic 
(bosonic) quarks, Sf (Sb) is the fermionic (bosonic) action, and P is a prefactor depending on 
all variables. Evaluating the partially quenched partition function ()35p at the saddle points, it 
may happen that the prefactor P is zero at some saddle points, such that these saddles must 
be disregarded. In this case the replica method can lead to wrong results. In appendix |Aj the 
occurrence of this phenomenon will be shown in the case of partially quenched theory of massless 
Wilson fermions. 

To summarize, we have computed the boundaries of the oscillating regions of the eigenvalue 
density of the baryon number Dirac operator within mean field chiral perturbation theory. As 
suggested in the example of section [II] we expect that the unphysical early onset of the average 
baryon number (see figure [2]) will disappear due to oscillations inside the regions marked by K on 
figure HI In the following section we show, using a RMM, that this is indeed the case. 
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IV. THE OSCILLATIONS ON THE MICROSCOPIC SCALE 



The strong oscillations have a period of order 1/V and we therefore need to compute the 
eigenvalue density on the microscopic scale in order to resolve them. Here we carry out this 
computation usinga random matrix model for the spectrum of the baryon number Dirac operator 



similar to that of 



24 



25[ . On the microscopic scale the random matrix model is fully equivalent to 



311 ] for a review of random matrix theory at nonzero chemical 



chiral perturbation theory, see e.g. 
potential. 

The RMM partition function with Nj quark flavors, all of mass m, in a sector of zero topological 
charge is given by 



7 N f ,n 
J N 



(W>; 



dfi{W) w{W) Yl det(D (m) + /x f ), 
f /=! 



(36) 



where [if are the chemical potentials of the different flavors (which at the end all will be set equal 
to fj,), If is a complex N x N matrix and w(W) is the Gaussian weight function 



(37) 



W (W) = exp[-NTrWW*]. 
The baryon number Dirac operator, Do(m), is given by 

iW ml at 

with Ijv the identity matrix. The equivalence between the RMM and xPT holds in the microscopic 
limit of QCD with nonzero chemical potential, where the microscopic quantities 



D (m) 



(38) 



m = mN and ft = (J,V~N 



(39) 



are kept fixed in the thermodynamic limit, N — > oo, see for example the review [3l|. After the 
standard Hubbard-Stratonovitvich transformation and saddle point integration one sees that the 
RMM is equivalent to the chiral partition function in microscopic limit, except for an overall factor 
of e N t N r ( see e.g. 



241]). under the identifications 



m = mN <-> mT,V and p 2 = [i 2 N <-> fi 2 F 2 V. 



(40) 



It turns out to be a challenging task to write the partition function (j36|) as an integral over the 
joint probability density function of the eigenvalues of I?o( m )- F° r this reason we will focus on the 
case of massless quarks, where the oscillations are expected to be maximally important and where 
we can immediately write down the joint probability density function of the eigenvalues of Dq. 
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A. The RMM for Massless Quarks 

The partition function (I36p simplifies in the massless limit, 

- N f 

Z* f ({n f };m = 0) = dn(W) w(W)Y[det(iW + n f )det(iW^ + fi f ). (41) 
J f=i 

For Nf = we immediately recognize this as the partition function for the ensemble originally 
solved by Ginibre 32|. If we denote the eigenvalues of W by {rjk}- The eigenvalues of the baryon 
number Dirac operator A are related to the eigenvalues rj of W as A = irj. The partition function 
can be written as 

N 



< /=o =A/- n/ d w^ |2 iA^)i 2 

(42) 

= AfT[ d 2 Vk e- N ^\ 2 det [K N ( Vi , V *)} 



k=l ' 

where A/o and M are normalization constants, A^(rj) is the Vandermonde determinant and the 
kernel Ajv(x,y) is defined from the orthogonal polynomials pk(x) as 

K N (x,y) = ^2pk(x)pk(y), with p k (x) = Ju~- ( 43 ) 



fc=o v 
The polynomials pk are orthonormal with respect to the Ginibre weight, i.e. 

d 2 r,e- N ^pMPj(v*) = (44) 



Methods to solve matrix models such as (I4ip . where characteristic polynomials are evaluated in 
a known ensemble, has been given in [33J. The trick is to use the identity 

K L 

A K+L (x) = A K (x)A L (y) Yl Y[( Xk ~ with y* = x k+i- (45) 

k=ie=i 

With this trick at hand it is straight forward to carry out the integration in (14ip . Up to irrelevant 
factors of normalization, we have 

N N J 

< / ({M/})=AT n / d ^ e- N ^\ 2 \A N ( V )\ 2 l{(i V k + ^vt + »f) 
k=i J f=i 

hf N r 

= a s° ^ 11 / dx k dy k A N+Nf (x)A N+Nf {y)e- Nx ^ 

Atf / ({»M) 2 £ = V (46) 



2 

k 

| N r 



Aiv>({*/i/}) 

~ A 77" V\2 ^ de ^r [KN+NfiWiWj)] 
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where the kernel K is defined in equation (|43p and (x, y) is defined as 



x k 



Vk 



for k < N 



for k > N 



and y k 



for k < N 
for k > N 



(47) 



From (146j) one sees that the eigenvalue density is given by 

N 



p N f(r l ,r l *,{ f , f }) = (Y / S 2 (Vk-r l ) 



k=i 



■ M'- 



-N\ n Y 



zg* (0*,}) A Nf ({i» f }y 



K N+Nf (i],i]*) 



KN+Nf (Vi WNf 



K N+Nf (in Np rf) ■■■ K N+Nf (ifi Nf ,ifi Nf ) 



(48) 



where J\f' is a normalization constant to be canceled by the normalization constant in Z^ f ({///})• 
The spectral density (I48p is written in terms of 77, the eigenvalues of W, but as mentioned previously 
these eigenvalues are related to the eigenvalues of the baryon number Dirac operator by A = ir/. 
The spectral densities for A and r] is related by a rotation by tt in the complex plane, such that 
baryon spectral density is easily obtained from (|48p . Some calculations are more compact in terms 
of the T) eigenvalues, and we will for this reason keep on writing expressions in terms of r] in the 
following section. 

The complex structure of the spectral density (|48p introduces oscillations with a microscopic 
period. In the following section we show that the oscillations are bounded within the region 
predicted in section ITO1 by mean field xPT. Futhermore we will explicitly demonstrate the crucial 
role of the oscillations. 



B. The average baryon density from the oscillating eigenvalue density 

For a chemical potential less than a third of the nucleon mass, the theory is dominated by pions 
and the partition function is therefore independent of the chemical potential at low temperature. 
In this section we will show how the ^-independence enters the RMM. Exactly as in the example 
of section UH the baryon number Dirac spectrum has a strong dependence on \i (cf. figure [5]), even 
though the average baryon density does not. As for the calculation for the chiral condensate 
the //-independence of the observable is obtained from a strongly oscillating region of the eigenvalue 
density. 

As already mentioned, the quenched spectral density (seen in the top panel of figure E]) is the 
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FIG. 5: Top panel: The quenched spectral density with N = 36. This density is identical to the Ginibre 
distribution. Bottom panel: The Nf = 1 spectral density with N = 36 and fi = 0.45. One sees how 
unquenching introduces a strongly oscillating region. The amplitude of the oscillations grows exponentially 
with N and the period is proportional to 1/N. The peaks of the oscillations has been cut; the maximum 
amplitude is more than a hundred orders of magnitude larger than the scale displayed. 



Ginibre distribution, 



The one flavor partition function f|46[) . is given by the kernel, 



(49) 



Z N f V) = #JV+l(*/W)> 



n 



(50) 



and agrees with the result derived in [25] by a different method. The one flavor spectral density is 
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given by ([45]). 

... ,9 r ff»r,ifn iiAjr*T , i (in n*\i 

(51) 



P J (77,7/ ,/i) = 2e |r " K N+1 (rj,rj) 



One can easily check that the spectral density (I5ip contain a circular region of radius /j with strong 
oscillation, exactly as predicted in section Hill The lower panel of figure [5] illustrates this. 

Here we will show that the oscillations due to the second term in (|51|) exactly cancel the 
unphysical early onset of the quenched baryon density. The baryon density is defined from the 
partition function as 

Vr, Nf=1 (iA - d ]ne\7 Nf=1 (u)] - dK N+l{¥,¥) I dp . , 

VUb W-d^ l0g[ZN (/i)J " K N+1 (ip,ip) • (52) 



We will now compute the average baryon number starting from the eigenvalue density, (|5ip using 

Vn N B !=1 (p)= [ d^-^—p^HviV*^), (53) 



and show that the oscillating nature is essential in order to get agreement with (|52p . 

After addition and subtraction of a kernel i£jv+l(v, V*) the one flavor spectral density given 
in (1511) can be written as 



\p i (rj,rj*,fj,)=w(ri,ri*)[K N+ i(ri,ri*)-K N+ i(ifj,,ri*)] 

- w(rj, rf) \ ^ N+1 V , {K N+l (ri, ip) - K N+ i(ip,ip))] , (54) 

with the weight w(r),rf) = e~ Nr,r] * . Integration over the first term of the spectral density, 

d 2 rj—^ [Kn+i(v,V*) ~ K N+ i(i(j,,ri*)] = -i / d 2 7]w(r], rf) V : Pfe(??*), (55) 

is zero due to orthogonality. In evaluation of the nonzero integral we will exploit that the eigenvalues 
of the baryon number Dirac operator come in pairs (ir], irf) combined with the fact that the weight 
w(x,y) and the kernel K n (x,y) only depend on the product xy to write, 

t r N;=i, , . f , 2 / *^K N+1 {ip,rf) K N+1 (rj,ip) - K N+1 (ip,,ip) 
Vn B } (p) = + i d rjw{rj,rj )— — — : 

.2 / *^K N+1 (r],-ifi) K N+ i(-iri,r]*) - K N+ i(ip,ip) 



i d 7]w(rj, rf)- 



K N+l (ip,ip) rj*+ip 
The first integral only has a contribution from r\ = ip and the second only from 77* = — i/j. The 
reason for this is that KN+i(ip,rj*) and Kn+i(t], —ip) are reproducing kernels in the space of 
polynomials of order less than N, and thus in a way are equivalent to delta functions 5 2 {rj — ip.) 
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and S (rj* + ifi), respectively. Using this we see that the average baryon number can be written as 



lim 



Tj— HfJ, 



lim^^ K N+1 (-ifi,7]*) - K N+ i(ifj,,in) 



+ 



<IKn + i (ifi, ifi) I dfi 



Tj* + ifi 



(57) 



This is precisely the result we found in (|52p . and the computation demonstrates directly the 
crucial role of the strong oscillations. Note that this result is true for all values of N. It therefore 
automatically applies in the microscopic limit, where the RMM is equivalent to xPT. 



n 



N f =l 
B 



/i 



fJ-c 

FIG. 6: Plot of the baryon density in the RMM for N = 36 and m = (which implies = 0) as a function 
of the chemical potential. The dotted line indicates the large TV limit. The baryon density is seen to be zero 
below the phase transition at fi c = 0.527 . . . 16] which is the RMM analog of a third of the nucleon mass. 



As mentioned above the RMM partition function and the xPT partition function differ by a 
trivial overall factor e N s^ 2 in the microscopic limit. If one want show that the /i-independence 
of average baryon density one need to multiply this trivial factor on the RMM partition function 
before further calculation. The one flavor baryon density is therefore given by 



.v, , d , rr ,N f =i, \ Nfl 2 dK N+ i(ifi,ifi)/dfi 



Vn B s (fi) = —log[Z N f (/x)e 



+ 27V/J. 



(58) 



To evaluate the average baryon density (|58p in the large N microscopic limit, we write the partition 
function as 



< /= V)e^ 



2 r(jv + i, -Nfi 2 



i 



r(7v + i) 



(JV-l)! 



dte N{lo g t~t)_ 



(59) 



19 



This integral can be evaluated with saddle point integration. There are two local maxima: One 



at t = 1 and one at t = —fj 2 . The critical value of the chemical potential ji C: which decides which 



maximum is the global one, is given by 



& + !og^c + 1 



0. 



(60) 



For < fi c the integral is dominated by t = 1 and partition function is independent of /j,, 
hence the average baryon number is zero. For \n\ > \i c the average baryon number is given by 
Vug 1 ~ 2N([i + This is shown on figure El The critical \i c is thus the RMM analog of 

a third of the nucleon mass [35|. Beyond the critical chemical potential the system is no longer 
dominated by the Goldstone bosons and the equivalence of the RMM and xPT breaks down. 

Finally, let us note that the spectrum of the baryon number Dirac operator has a 27rnT- 
periodicity along the imaginary axis, where T is the temperature. This is a manifestation of 
the fact that if ip is an eigenfunction of .Do(m) with eigenvalue A, then e %UinX °ip is an eigenfunction 



10j. Promt 



with eigenvalue A + iu: 
uj n = 27rmT, where n is an integer 
chemical potential, see 



le boundary condition of ip in the time direction it follows that 



36]. This periodicity is related to the periodicity of imaginary 



37], but is not seen in the microscopic limit, since the eigenvalues scale as 



A ~ T 2 



V. CONCLUSIONS 



In QCD with nonzero chemical potential the sign problem manifest itself in the physical ob- 
servables. This is true for the average baryon number as well as for the chiral condensate. If one 
neglects the phase factor of the fermion determinant the baryon density will have an unphysical 
early onset at fi = m n /2 and the chiral condensate will rotate into a pion condenssate. 

In this paper the correct physical behavior of the unquenched average baryon number has 
been linked to strong oscillations of the eigenvalue density of the baryon number Dirac operator. 
The oscillations have a period on the microscopic scale and were shown to be dominant within 
bounded regions of the complex eigenvalue plane. The mechanism which links the oscillations 
of the baryon number Dirac spectrum to the baryon number density is in exact correspondence 
with the mechanism between the chiral condensate and the Dirac spectrum. This shows the general 
nature of this mechanism, when a sign problem is present. The boundaries of the oscillating regions 
of the baryon number Dirac spectrum where computed within mean field chiral perturbation theory 
and using a random matrix model it was shown exactly how the oscillations are responsible for the 
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//-independence of the average baryon density. 

In conclusion, we have shown that the mechanism which links physical observables in un- 
quenched QCD to strong oscillations of the corresponding eigenvalue density is not restricted 
to the chiral condensate; it also holds for the average baryon number. This solves the Sliver Blaze 
problem [10] for the baryon number Dirac operator. The oscillations have a period of order 1/V and 
an amplitude that grows exponentially with V. Moreover, in both cases one must integrate over 
at least V periods of the oscillations to approach the correct physical behavior. This demonstrates 
how severe the sign problem is for \x > m n /2, see also 

It is an appealing mathematical challenge to generalize the computation carried out here within 
the random matrix framework to the case of a nonzero quark mass. The direct supersymmetric 
computation of the microscopic spectral density is extremely demanding in this case, but perhaps 
it can be simplified through the use of integrable structures as in {39!. l^ol]. 

It would be interesting to extend this study to the ldQCD baryon number Dirac operator. 



Finally, it wou 
approach of 



411. 



d also be most interesting to study the results of the present paper within the 
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Appendix A: Partially Quenched Wilson Fermions 

In section imi it was mentioned that saddle point integration within the replica method may give 
saddle points that should be neglected. This happens when the corresponding saddle point in the 
partially quenced theory has a prefactor identical to zero. In this appendix we show an example 
of this phenomenon in the effective theory of QCD with Wilson fermions. 

The partition function describing QCD with Nf flavors of Wilson fermions in a sector with zero 
topological charge is 

Z N JM,a)= f dUe~ s(U ' M ' a \ (Al) 

Ju(N f ) 

In the microscopic limit the action is given as 

S(U, M,a) = -^VTt M(U + W) 

- a 2 W 6 V[Ti(U + U ] )] 2 - a 2 W 7 V[Tv(U - U ] )] 2 - a 2 W 8 V Tv(U 2 + C/ t2 ), (A2) 
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where M is the mass matrix, a is the lattice spacing, S is the chiral condensate and W% are low 
energy constants which determine the leading order discretization errors of Wilson fermions 42j. 



Here we will focus on the special case where W% = Wj = 0. The scaled variables, M = MT,V and 
a| = a 2 WsV, are kept at order unity in the microscopic limit. One finds the partially quenched 
partition function 4314451] 

Z Nf+lll (M,a 8 ) = j dUe ^ti(u-u-i)+alSHU*+u-*)_ (A3) 
The integration is over the maximum Riemannian submanifold of Gl(Nj + 1|1). For Nj = we 



can parameterize the graded manifold as 



46] 



(e ie 0) 




( a \ 




exp 




\ e s J 




\P 0/ 



U=\ exp , (A4) 

\ e s ) \P Oj 

where a and /3 are Grassmann variables. Setting M = ml and integrating out the Grassmann 
variables gives 

-i poo rn 

Z m (m,a 8 ) = — ds d0P(m,a,9,s)e- s f ( - m ' a8 ' 0) e- Sb< - m ^' s \ (A5) 

27r J-OC J~TT 



with 



Sf(m, a 8 , 6) = + msin# + 2ag cos 26, (A6) 

Sb(m, as, 9) = + im sinh s — 2a\ cosh 2s, (A7) 

,-./ r\ \ rh . n ,rh , . 
F{m, a 8 , t),s) = — — smd + i— smns 

+ 2a§ (cos 26 + cosh 2s + 2 cos 6 cosh s — 2i sin 6 sinh s) . (A8) 



Numerically one can easily verify that the partition function (1A5|) equals one as it should, in- 
dependently of rh and as, since we work at equal fermion and boson sources. Here we want to 
calculate (|A5P with saddle point integration. The saddle points are given by 

s = 0, cos 6 = (A9a) 

* 2 
Tfi 

s = 0, sine = ( A9b ) 

8a 8 

In general both saddle points will contribute to the partition function, but note that the prefactor 
P evaluated at the saddle point (|A9aP is 

771 

P(m,as,6 = -^,s = Q) = -, (A10) 
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such that the prefactor becomes zero in the massless limit. Calculating the partition function (|A5 
in massless limit with saddle point integration one only needs to include the saddle point (|A9bj) . 

P(a,9,s)e- s t( a S'^e- s »( as ' s *> 



1. 



(All) 



=s=0 



\S'}(a 8 ,e)S^a 8 ,s)\ 

This behavior is in complete analogue to the behavior of the massless theory mentioned in sec- 
tion [Till In such cases the mean field replica method can lead to wrong results. 
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